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Abstract. In the light of the "no-hair" conjecture, we revisit stable black holes in 
su(iV) Einstein- Yang-Mills theory with a negative cosmological constant A. These 
black holes are endowed with copious amounts of gauge field hair, and we address the 
question of whether these black holes can be uniquely characterized by their mass and 
a set of global non-Abelian charges defined far from the black hole. For the su(3) case, 
we present numerical evidence that stable black hole configurations are fixed by their 
mass and two non-Abelian charges. For general N, we argue that the mass and TV — 1 
non-Abelian charges are sufficient to characterize large stable black holes, in keeping 
with the spirit of the "no-hair" conjecture, at least in the limit of very large |A| and 
for a large subspace of the total space of stable black holes. 



1. Introduction 

According to the "no-hair" conjecture p], equilibrium black holes are extraordinarily 
simple objects, characterized completely by their mass, angular momentum and charge. 
These three quantities, as well as being parameters in the Kerr-Newman metric, are 
also physical, global charges, which could, at least in principle, be measured far from 
the black hole. For stationary, asymptotically flat black holes in four- dimensional 
Einstein-Maxwell theory, the "no-hair" conjecture has been proved (see for example 
[2J for reviews). It is perhaps unsurprising that, if one of the above assumptions 
(asymptotically flat space-time, four space-time dimensions, Einstein-Maxwell theory) is 
relaxed, then black hole uniqueness no longer holds. Here we are particularly interested 
in matter theories other than Einstein-Maxwell, and we focus on the four-dimensional 
Einstein- Yang-Mills (EYM) model, on which there is now an extensive literature (see, 
for example, [3111] for reviews). 

The original discovery of soliton [5] and black hole [6j solutions of su(2) EYM 
theory in four- dimensional, asymptotically fiat space, seemed, at least at first, to provide 
counter-examples to the "no-hair" conjecture. The black hole solutions have no global 
charges and are indistinguishable from the usual Schwarzschild black holes at infinity. 
However, both the black hole and soliton solutions are unstable [7] , so while the "letter" 
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of the "no-hair" conjecture (which, in its original form, says nothing about stabihty) is 
violated, its "spirit" remains intact, because stable black holes do seem to be uniquely 
characterized by global charges. These results led Bizon to reformulate the "no-hair" 
conjecture as follows [8]: 

Within a given matter model, a stable stationary black hole is uniquely 
determined by global charges. 

For four- dimensional EYM theory in asymptotically flat space, this result holds at least 
for spherically symmetric black holes with a purely magnetic gauge field, as the unique 
stable solution in this case is the Schwarzschild black holqj. For asymptotically fiat, four- 
dimensional EYM black holes with both an electric and magnetic field, the situation is 
less clear, with the recent discovery [10] of stable dyonic black holes with a non- vanishing 
electric charge but zero magnetic charge. 

When a negative cosmological constant A is introduced into the model, so that 
the space-time is asymptotically anti-de Sitter (adS) rather than asymptotically fiat, 
stable soliton [TT] and black hole solutions [12] of 5u(2) EYM theory exist. For 5u(2) 
gauge group, the purely magnetic gauge field is described by a single function u. Purely 
magnetic solutions have been proven to be stable [HI [121 13] 5 when the function u 
has no zeros, provided that |A| is sufficiently large. Unlike their asymptotically flat 
counterparts, these solutions do have a global magnetic charge. 

Enlarging the gauge group to 5u{N), the purely magnetic gauge field is now 
described by — 1 functions Uj. For any fixed A^, the existence of stable, purely 
magnetic, soliton and black hole solutions (for which all the uj have no zeros) has been 
proven provided that |A| is sufficiently large [HI [151 UHl [II]- These stable black holes 
have A^ — 1 independent gauge degrees of freedom and it turns out (see [I6] and section 
Hj) that they are described by A^ + 1 independent parameters, for gauge group su{N). 
We conclude [H] that there is no limit to the amount of stable gauge field hair with 
which a black hole in adS can be endowed. 

Our purpose in this paper is to revisit these su{N) EYM black holes in the light of 
the "no-hair" conjecture as reformulated by Bizon. In particular, we investigate whether 
these black holes are uniquely characterized by global charges. After a brief review of 
the su{N) EYM model and its black hole solutions in section [2], we proceed, in section [3l 
to construct global charges for the 5u{N) gauge field. We follow a two-pronged approach 
to examine whether these charges uniquely characterize the black hole solutions, firstly, 
in section HI performing numerical investigations of the solution space, and then, in 
section [6], giving an analytic argument that the black holes are uniquely characterized 
by global charges, at least for sufficiently large |A|. Our conclusions on the consequences 
of this result for the "no-hair" conjecture are presented in section [71 

I It has been proven that aU sphericaUy symmetric, four-dimensional, black hole solutions of EYM 
theory in asymptotically flat space with a purely magnetic gauge field are unstable, for all gauge 
groups [9]. 
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2. Hairy black holes in 5u{N) Einstein- Yang-Mills theory in anti-de Sitter 
space 

In this section we outline the salient features of su(A^) EYM theory in anti-de Sitter 
space, and the black hole solutions found numerically in [I5j, whose existence, for 
sufficiently large |A|, was proven in [16]. We also discuss the thermodynamic properties 
of the black holes, and use the boundary counter-term formalism [18] to compute their 
mass. 



T,. = F^,F:' - -g,,FlF^"^, (3) 



2.1. Ansatz and field equations 

We consider four- dimensional 5u{N) EYM theory with a negative cosmological constant, 
described by the following action, in suitable units: 

1 ^4. 



^EYM = -J d^x^ [R-2A- Tr F^^F^"] , (1) 

where R is the Ricci scalar, A the cosmological constant and Tr denotes a Lie algebra 
trace. Throughout this paper, the metric has signature (— , +, +, +) and we use units 
in which AirG = 1 = c. In addition, we fix the gauge coupling constant g = 1. In this 
article we focus on a negative cosmological constant, A < 0. Varying the action ([T|) 
gives the field equations 

=D^F/ = V^F/+[A^,F/]; (2) 

where the YM stress-energy tensor is 

1 

4' 

with summation over the Lie-algebra index a understood, so that T^i, involves a Lie- 
algebra trace. The Yang-Mills gauge field F^y is given in terms of the gauge potential 
by 

F^,u = d^Ay - d^A^ + [A^, Ay] . (4) 

In this paper we are interested in static, spherically symmetric black hole solutions 
of the field equations ([2]), and we write the metric in standard Schwarzschild-like co- 
ordinates as: 

ds^ = -fiS^ dt^ + /i"^ dr"^ + dO'^ + sin^ d(j)'^, (5) 

where the metric functions /i and S depend on the radial co-ordinate r only. In the 
presence of a negative cosmological constant A < 0, it is convenient to write the metric 
function /i as 

2m(r) Ar2 

^ r = 1 — . (6) 

r 6 

We emphasize that in this paper we are considering only spherically symmetric black 
holes and not topological black holes which have been found in the su(2) case 
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With a suitable choice of gauge, we take the purely magnetic 5u{N) gauge potential 
to have the form [20| 



A = ^{C -C^) de [{C + C^) sine + D cose] d(f), (7) 

where C and D are (A^ x A^) matrices and is the Hermitian conjugate of C. The 
constant matrix D takes the form: 

D = Diag (A^ - 1, - 3, . . . , -AT + 3, -A^ + 1) , (8) 

and the matrix C is upper-triangular, with non-zero entries only immediately above the 
diagonal: 

Cj,j+i = ujj{r), (9) 

for j = 1, . . . , A^ — 1. The gauge field is therefore described by the A^ — 1 functions ujj{r). 
The derivation of the ansatz ([7]) uses the Yang-Mills equations and assumes that all the 
u)j{r) are not identically zero (see, for example, [21] for the possibilities in asymptotically 
flat space if this assumption does not hold). We comment that our ansatz (I7j) is by no 
means the only possible choice in su{N) EYM. Techniques for finding all spherically 
symmetric 5u( A^) gauge potentials can be found in [22] , where all irreducible models are 
explicitly listed for A^ < 6. 

With the ansatz ([7]), there are A^ — 1 non-trivial Yang-Mills equations for the A^ — 1 
gauge field functions uj: 

( 2kr^\ 
r^jjooj + f 2m - 2r^pe — j ^ + WjUj = 

for j = 1, . . . , A^ — 1, where a prime ' denotes d/dr, 

1 ^ 

=^.Y.[i^^--^U-N-i + 2jY], (11) 

W^. = l-^J + ^K-i + ^J+i), (12) 
and = ujn = 0. The Einstein equations take the form 

5" 2G 

m' = fiG + r^Pe, — = — , (13) 

b r 

where 

N~l 

G=Y,uf. (14) 

The field equations (ITU| [T^ are invariant under the transformation 

uj.j{r) -ujj{r) (15) 

for each j independently, and also under the substitution: 

j^N-j. (16) 



(10) 
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2.2. Boundary conditions 

Our primary interest in this paper is black hole solutions of the field equations ( fTO| 
[T3i) . However, we will need to consider solitons in section I6l2] to ensure that solitons 
cannot be mistaken for black holes by measuring global charges at infinity. The field 
equations are singular at the origin, at an event horizon r = rh ii there is one, and at 
infinity r — t- oo. Boundary conditions therefore have to be specified in a neighbourhood 
of these singular points. Local existence of solutions of the field equations satisfying the 
boundary conditions outlined below is proven in fl6j . 

2.2.1. Origin The boundary conditions at the origin are more complicated than near 
the event horizon or at infinity. The full form is derived in detail in pTO] (following the 
analysis of [23] for the asymptotically fiat case). Here we simply state the basic features 
which are needed for our analysis in section 16.21 Near the origin, the field variables have 
the following form: 

m(r) = m-^r^ + O(r^); 
Sir) =SQ + S2r^ + 0{r^)- 

uj,{r)= ±[j(iV-j)]^ + 0(r2). (17) 

To fully specify the form of the gauge field in a neighbourhood of the origin, a 
complicated power series has to be developed, up to O(r^), the details of which can be 
found in [16j but which are not necessary for our purposes in this paper. 

2.2.2. Event horizon For black hole solutions, we assume that there is a regular, non- 
extremal event horizon at r = r/j, where /i(r) has a single zero. This fixes the value of 
m(r/i) to be 

m(r.) = y - ^- (18) 

We assume that the field variables a;j(r), m(r) and S{r) have regular Taylor series 
expansions about r = r^. 

m{r) = m{r-h) + m'{rh) (r — r^) + O [r — r^f ; 

ujj{r) = Uj{rh) + u;j{rh) (r - r/,) + O (r - rtf ; 

S{r) = S{rh) + 5'(r,) (r - r,) + O (r - r,) . (19) 

Setting /i(rft) = in the Yang-Mills equations ( ITOl) fixes the derivatives of the gauge 
field functions at the horizon: 

Wj{rh)uj{rh) 



^■{rn) = "nnr-,,nj ^^^^ _ ^20) 

2m{rh) - 2rfj)g{rh) 3-^ 

Therefore the expansions f lT9|) are determined by the + 1 quantities Uj{rh), rt, S{rh) 
for fixed cosmological constant A. For the event horizon to be non-extremal, it must be 
the case that 

2m'(rO = 2rlpg{rh) < 1 - Arl (21) 
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which weakly constrains the possible values of the gauge field functions ujj{rh) at the 
event horizon. Since the field equations f lTOl [T3|l are invariant under the transformation 
f lTSl) . we may consider ojjirh) > without loss of generality. 



2.2.3. Infinity At infinity, we require that the metric ([5]) approaches adS, and therefore 
the field variables ujj{r), m{r) and S{r) converge to constant values as r — )■ oo. We 
assume that the field variables have regular Taylor series expansions in near infinity: 

m{r) = M + O [r-^) ] 5(r) = 1 + O (r^^) ; 

U0j{r) = Uj,oo + c,r~^ + O {r'^) . (22) 

We have included the 0{r~^) terms in a;j(r) as they are central to our analysis in 
section HI If the space-time is asymptotically fiat, with A = 0, then the values of ojj^oo 
are constrained to be 



.,,00 - ±^j{N-j). (23) 

This condition means that the asymptotically flat black holes have no magnetic charge 
at infinity (see section [3] for the definition of magnetic charges). Therefore, at infinity, 
they are indistinguishable from Schwarzschild black holes. However, if the cosmological 
constant is negative, then there are no a priori constraints on the values of Uj^oo- In 
general, therefore, the adS black holes will be magnetically charged. In section [3] we will 
construct appropriate non-Abelian charges. 

2.3. Embedded solutions 



The field equations (jlOl fT3i) are non-linear and coupled, but they do have two analytic, 
trivial solutions. 

Schwarzschild-adS Setting 



uj,{r) ^ ±^j{N-j) (24) 

for all j gives the Schwarzschild-adS black hole with 

m{r) = M = constant. (25) 

Reissner-Nordstrom-adS Setting 

LOjir) = (26) 

for all j gives the Reissner-Nordstrom-adS black hole with metric function 

2M Q2 ^2 

Mr) = l- — + ^-— , (27) 

where the magnetic charge Q (see section [3] for a definition of this quantity) is fixed 
by 

Q^ = lN{N + l)iN-l). (28) 

Only for this value of the magnetic charge is the Reissner-Nordstrom-adS black hole 
a solution of the field equations. 
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As well as these effectively Abelian embedded solutions, there is also a class of 
embedded su(2) non- Abelian solutions, given by writing the — 1 gauge field functions 
u}j{r) in terms of a single function uj{r) as follows: 



It is shown in [16] that, with a suitable rescaling of the other field variables, the field 
equations ( jTOl [T3|) reduce to the su(2) field equations for the function co. Therefore any 
su(2), asymptotically adS, EYM black hole solution can be embedded into su{N) EYM 
to give an asymptotically adS black hole. 

2.4- Properties of the su{N) solutions 

The properties of soliton and black hole solutions of the field equations ( fTOl [T3l) have 
already been studied in detail elsewhere [H HH [151 HH], therefore here we simply 
summarize the salient features required for our subsequent analysis. 

For large |A|, numerical investigations [HI |T2] find both soliton and black hole 
solutions for which all the gauge field functions 0Jj{r) have no zeros. It has been proven 
[16] that, for fixed rh and Uj{rh), then black hole solutions for which all the O0j{r) have 
no zeros exist for all sufficiently large |A|. In view of these results, the focus in numerical 
work [15] has been on properties of the phase space of solutions for fixed r^ and different 
values of A. In this work we are interested in the properties of the phase space for fixed 
A and varying horizon radius r^. 

We find that, for sufficiently large |A|, all numerical solutions for varying are 
such that the gauge field functions tjj{r) have no zeros. For example, in figure [1] we 
show the phase space of solutions for su(2) black holes with A = —10 and varying Vh- 
Corresponding phase space studies for larger A^ with A fixed and rh variable are more 
complex because of the number of parameters involved, but in each case investigated we 
find similar results. Phase space plots for N > 2 can only be produced by fixing some 
of the parameters, see, for example, figure [2l 

At least some of these black hole solutions for which all uj have no zeros are stable 
[m [17]. For su(2), black hole solutions for which the single gauge field function u{r) 
has no zeros are shown to be stable if w(r) > 1/ -\/3 everywhere on and outside the event 
horizon [12]. For N > 2, stability under linear, spherically symmetric perturbations 
can be proven for sufficiently large |A|, for su{N) solutions satisfying the inequalities 



and which are sufficiently close to stable embedded su(2) solutions. The inequalities (1501) 
have to be satisfied for all values of r > rh- In figure [2] we plot the region of phase space 
in the 5u(3) case, with A = —3 and rh = 1, for which the inequalities fl30|) are satisfied 
on the black hole event horizon and for which both Ui{r) and O02{r) have no zeros. For 
at least some of these solutions, the inequalities ( 130|) are satisfied for all r >rh [13| [T7] . 
However, the inequalities fl3Up are satisfied only in a comparatively small region of the 



;,(r) = ±^j{N-j)u{r) Vj = 1, . . . , AT - 1. 



(29) 



[lain]: 




(30) 
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Figure 1. Phase space of su(2) solutions with A = — 10 and varying r/j. The insert is 
a blow-up for smaller values of r^, so that the structure can be seen. All the numerical 
solutions we find are such that U!{r) has no zeros. The light grey region corresponds to 
values of a;(r;i) such that the condition ()21|) is satisfied, but we do not find a numerical 
solution. Inside the light grey region, in the region labelled "n = 0" , we find numerical 
solutions for which uj{r) has no zeros. Outside the light grey region, the condition (|2T|) 
for a regular event horizon is not satisfied. 




Figure 2. Phase space of su(3) solutions with A = —3 and — 1, with ni the number 
of zeros of the gauge field function uji{r) and n2 the number of zeros of L02{r)- We 
plot only the region in which the inequalities pop are satisfied on the event horizon. 
The light grey region indicates solutions for which both the gauge field functions uJi{r) 
and (^2(7') have no zeros. The dark grey region corresponds to values of ijJi{rh) and 
W2(r'?i) for which the inequalities pOI) are satisfied on the event horizon, together with 
the condition (pT|) for a regular event horizon, but for which we do not find solutions. 



phase space for which there are nodeless solutions (compare the size of the region in 
figure |2] with figure 2 in [16], where the entire phase space for 5u(3) black holes with 
A = — 3 and r/j = 1 is plotted). Furthermore, due to the algebraic complexity of the 
gravitational sector perturbation equations, stability for this sector can only be proven 
for sufficiently large |A| and for su(iV) solutions sufficiently close to stable embedded 
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Figure 3. Entropy-temperature curves for su(2) black holes with A = —10. For each 
curve, the values of the gauge field function uj(r) at infinity is fixed. The curve with 
a;(oo) = 1 corresponds to embedded Schwarzschild-adS black holes. 



su(2) solutions, with, as is usually the case in this kind of proof, little indication of 
how large "sufficiently large" is, or how close "sufficiently close" is [17]. Our focus in 
this article is the characterization of stable su(iV) black holes. We therefore restrict 
attention to those black holes for which the inequalities fl30|) are satisfied, at least on 
the event horizon. It is likely that our results include some black holes which are in fact 
unstable. 

2.5. Thermodynamics of su{N) EYM black holes 

While the "no-hair" conjecture as outlined in the introduction is primarily concerned 
with the characterization of black holes which are c/assica//?/ stable, it nonetheless makes 
sense to consider the thermodynamic stability of the black holes. The entropy of su{N) 
EYM black holes is given, as usual, by one-quarter of the area of the event horizon, 
while their Hawking temperature Th is 

Th = -^{1- 2m'(r,) - Ar^) 5(r,), (31) 

where we emphasize that S{rh) is a metric function ([5]), and not the entropy of the black 
hole. For su(2) black holes, the thermodynamics has already been studied for A = —3 
in [21]. Their results are similar to ours in figure [3] for the A = — 10 case. In figure [3l 
we plot entropy as a function of Hawking temperature for fixed ci;(oo) (we will see in 
section [3] that this corresponds to fixing the magnetic charge of the black hole). For each 
value of u;(oo), there are two branches of black hole solutions. Firstly, there is a branch 
with small entropy and negative specific heat, so that the black holes on this branch 
are thermodynamically unstable. Secondly, there is an upper branch of black holes with 
larger entropy and positive specific heat, corresponding to thermodynamically stable 
black holes. This is the behaviour found by [23] when A = —3, and we found similar 
behaviour for other values of A. The curve in figure [3] with u{oo) = 1 corresponds to 
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Figure 4. Entropy-temperature curves for su(3) black holes with A = —3. For each 
curve, the values of the gauge field functions Wi(r) and uj2{r) at infinity are fixed. 
Curves with a;i(oo) = a;2(oo) correspond to embedded su(2) solutions. 




embedded Schwarzschild-adS black holes. As predicted in [25] for generic hairy black 

holes in asymptotically fiat space, we see from figure |3] that the non-Abelian su(2) black 
holes have lower temperatures than the embedded Schwarzschild-adS black hole with 
the same entropy (and hence horizon area). Similar behaviour is observed for su(3) 
black holes, as can be seen in figure |H In figure 4, again anticipating the results of 
section [31 we fix the values of the two gauge field functions Co'i(r) and U2{r) at infinity, 
and then plot the curve of entropy as a function of temperature. If (X'i(oo) = ^2(00), 
then we have embedded su(2) solutions and the curves have a very similar shape to 
those in figure [3l with two branches of solutions. In this case we are interested in the 
upper branch of solutions which are thermodynamically stable. For ui{oo) 7^ ci;2(oo), 
we have genuinely su(3) solutions, and in this case the curves appear to have just one 
branch of solutions, which have positive specific heat and so are thermodynamically 
stable. We are not able to say conclusively whether the fact that we have been unable 
to find a thermodynamically unstable branch of solutions is due to numerical difficulties 
or whether such a lower branch of solutions does not in fact exist. 

2.6. Mass ofsu{N) EYM black holes 

The mass of the non-Abelian black holes is readily computed using the background 
counter-term formalism [18]. As is well known, the gravitational action for 
asymptotically adS space-times typically diverges as the boundary dAi of the region 
(over which the Lagrangian is integrated) is taken to infinity. This divergence is removed 
by the addition of boundary counter-terms to the gravitational action, which do not alter 
the bulk equations of motion. As we are working in four space-time dimensions, the 
boundary counter-terms which are sufficient to yield a finite bulk gravitational action 
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are 



i-ct 



d'xV^i^ + in], (32) 



dM 



2 



where i = y ~f is the adS length and the boundary metric 7 has Ricci scalar TZ. The 
boundary stress tensor, resulting from the variation of the total gravitational action 
'S'grav (which is the usual gravitational action plus the counter-terms f l32|) ) is then [18] 



^5 = = I f^M. - ©7.. - ^7.. - iG,. ] , (33) 



where 9^;^ is the extrinsic curvature of the boundary (with trace 9), and is the 
Einstein tensor of the boundary metric. Only the component T^f is required to compute 
the mass of the black holes. To leading order, we find that 

TS = f (34) 

where M = limr^oo '^('") (I22|) . so that the mass of the solutions is simply 

frT^f d^x = AttM. (35) 



B 
tt 

dM 

From here on we will therefore use the variable M to denote the "mass" of the black 
holes. Our results here are in complete agreement with those in [21] , where it was found, 
for the su(2) case, that the non-Abelian gauge fields only contribute to the boundary 
stress tensor at next-to-leading order, and therefore do not contribute directly to the 
mass of the solutions. Of course, they do contribute indirectly to the mass as they affect 
the value of M through the Einstein equations. 



3. Defining charges for 5u{N) EYM black holes 

We now turn to the definition of global charges for su(A^) EYM black holes. Since 
the rank of the su(iV) Lie algebra is — 1, we expect these black holes to carry — 1 
conserved charges, which, at least in principle, could be measured at infinity, far from the 
black hole. The definition of charges for general non-Abelian gauge fields is non-trivial 
because of the need for observable charges to be gauge-invariant, coupled with the fact 
that the gauge field F^i, is not itself gauge-invariant in the non-Abelian case (unlike the 
situation for Abelian gauge fields). Methods for defining global, gauge-invariant, charges 
for non-Abelian gauge fields have been devised by a number of authors [26], ET] EH] 129] . 
While the quantization of non-Abelian charge has been studied [27], here we take a 
purely classical approach. 

For an 5u(2) gauge field, the single magnetic charge was defined in [30] as 

Q = ^[ \/W^dedcP, (36) 

where, under the square root, we have a Lie algebra trace over the Fg^f, component of 
the gauge field, corresponding to a sum over the Lie algebra index a. The integral is 
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taken over the sphere at infinity. The disadvantage of the formula fl36p is that there is 
no natural generalization yielding A^ — 1 charges in the su{N) case. An alternative su(2) 
charge of the form 

Q = l-iol (37) 

has been considered by a number of authors (see, for example, [21]), and has the 
advantage of readily extending to the larger gauge group. In particular, in this section 
we use the definitions of [28] (similar quantities were also defined in [29]) to construct 
magnetic charges. The approach of [28] has the advantage of yielding expressions which 
can be easily applied to our solutions. 

Following [28], we define gauge invariant magnetic charges as follows. Let X be an 
element in 21, the Cartan sub-algebra of the 5u{N) Lie algebra. Then, for each X, a 
charge Q{X) is defined by: 

Q{X) = ^ sup A; (x, I g-^Fg\ . (38) 

Here the integral is taken over the sphere at infinity, and the supremum is taken over 
all group elements g{r) of the form 

g{r) = exp[f{r)^ (39) 

where /(r) is a scalar function of r and E is a constant element in the Lie algebra. 
In (1381) . we have used k{X,Y) = TrjadXadF} which is the Killing form, with adX 
denoting the adjoint representation of the Lie algebra element X. It is shown in [28] that 
the integrand in fl38|) takes its maximal value when g~^Fg is in the Cartan sub-algebra 
21, and therefore we may restrict attention to those g{r) for which this is the case. 

Since we are integrating over a sphere at infinity, the integral in (l38ll simplifies to 

Q{X) = ^k (x, I g-'Fe^g dO dA . (40) 

From the gauge potential ansatz ([7]), the required component of the field strength is 
found to be 

Fe^ = -^{[C,C'']-D)sme, (41) 

where D is the constant matrix given in ([H]). From the form of the matrix C OH]), it is 
straightforward to show that 

[C,C^] =diag{ujf,uj^,-ujlujl-ul...,-uj%_,}. (42) 

We find that Fg^ (HTi) is an element of the Cartan sub-algebra 21. To see this, we begin 
by defining the following generators of the Cartan sub-algebra 



Hk = --^_==diag <^ M^l, -fc, 0, 0, . . . , . (43) 

k entries 
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Then Fg,^ can be written in terms of the generators as follows: 

F,,^^2^Y.V2kWTT)(f-^-l)H.. (44) 



, k k+1 
k=i \ ' / 

Since Fq^ is itself an element of the Cartan sub-algebra 21, we can set g{r) = e, the 
identity element, in fHOj) . For each element of the Cartan sub-algebra X G 2t, equation 
( l40l) yields a charge Q{X), however only N — 1 of these will be independent because the 
rank of the su{N) Lie algebra is N — 1. We therefore need a consistent way of selecting 
an appropriate basis of charges Qi, ioi i = 1, . . . , N — 1. 

First of all, we require our su{N) charges (140!) to reduce to (multiples of) the su(2) 
charge (|37l) for embedded su(2) solutions. Secondly, we fix an overall normalization 
constant by defining a total "effective" charge Q (reminiscent of the alternative 5u(2) 
charge f l36|) ). as follows. As r — )■ cxo, the metric function m(r) takes the form 

m{r) = M - — + 0{r-^), (45) 
with Q given by 

Af-l ^ N 

Q' = Y.Q'=2^ Koo - -N-l + 2j)\ (46) 

i=i j=i 

Therefore Q plays the same role in the metric as the usual Abelian charges. For 
embedded Reissner-Nordstrom-adS black holes (126|) . the expression ( H6l) reduces to 

Q2 = liv (iV - 1) (AT + 1) . (47) 
6 



With these two constraints, the natural basis of charges for su{N) is 

(4J 



_ y/j (j + 1) 
— 



2 2 

1 _ _^ '^i+i.oo 



J J + 1 

for j = 1, . . . , X — 1. For the su(2) case, there is just one charge Qi = 1 — (X'(oo)^, 
as required. For the 5u{N) case and embedded 5u(2) black holes with Uj{r) = 
- j)u{r), the formula (HE]) reduces to 



Q, = y^I^^[l-oo{oor] (49) 

so that all the su{N) charges are proportional to the su(2) charge, again as required. It 
is also straightforward to check that fj46|) holds. 

From (jUD, it is clear that the charges are determined uniquely by the values of the 
gauge field functions Uj at infinity. What is not so immediately apparent is that the 
converse is also true: the values of the gauge field functions at infinity can be determined 
(up to an overall irrelevant sign) from the charges Qj fHHj) . In particular, from fj48l) . it 
can be shown that 

N-l 

^loo =J{N- J) - jV2j2 Ji—r. - (50) 

k= 



- v/fc(^Tiy 
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4. Characterizing 5u{N) EYM black holes for large |A| - numerical work 

Having defined a set of global charges for our 5u{N) EYM black holes, we now address 
the question of whether these charges, together with the mass M defined in section 12.61 
and the negative cosmological constant A, are sufficient to characterize stable 5u{N) 
EYM black holes. We begin, in this section, with a numerical investigation, before 
turning, in section [6l to analytic arguments. 

In [16], two local existence theorems are proved, for solutions in a neighbourhood of 
the event horizon and infinity. Near the event horizon, it is proven that the parameters 
A, rh and ujj{rh) ( |T9|) determine the solutions in a neighbourhood of the event horizon, 
and, furthermore, that these local solutions are analytic in these parameters and the 
variable r. Therefore iV + 1 parameters are required to completely specify the black 
hole solution^. 

Near infinity, the situation is different. It is proven in [16] that 2N parameters are 
required to determine the local solutions in a neighbourhood of infinity, namely A, M, 
u)j^oo and Cj (122|) . Since we know that the black hole solutions form an + 1 parameter 
family, it is clear that only + 1 of these 2N parameters are independent, but the 
question is, which A^ + 1 parameters can we take to be a basis? Ideally we would like 
the A^ — 1 charges Qj fj48|) . defined in the previous section, together with the mass M 
and cosmological constant A, to be a suitable basis. From the analysis in the previous 
section, the charges Qj are completely specified by the values of the gauge field functions 
Uj^oo at infinity, so equivalently we would like the constants cj in (|22l) to be determined 
by A, M and Uj^oo- 

We begin our numerical investigations with the simplest case, namely su(2) black 
holes. In figure [5] we plot ci against the mass M and charge Qi for stable su(2) black 
hole solutions when A = — 10 (that is, we consider only solutions for which the gauge 
field function u{r) has no zeros and u{rh),uJoo > !)• It is clear from the surface in figure 
[5] that Ci is a single- valued function of M and Qi- We find similar results for other 
values of A. For further evidence that M and Qi therefore characterize 5u(2) black 
holes, in figure [6] we plot M as a function of Qi for various r^, with A = — 10. We have 
found no evidence that the {Qi, M) curves for different r^ ever intersect. Therefore, for 
each value of [Qi, M) for which there is a black hole solution, that solution is uniquely 
characterized by Qi and M. 

For su(3) and larger gauge groups, it is more difficult to see graphically whether 
the Cj are uniquely determined by M and the charges Qj, simply because there are more 
variables to plot. However, all the numerical evidence we have been able to gather in 
the su(3) case does indicate that the Cj are uniquely determined by M and Qj. To 
illustrate this, we consider 5u(3) black holes with A = —3, such that the inequalities 

§ Strictly speaking, there is one further parameter needed near the horizon, namely S{rh). However 
this parameter is fixed by the requirement that S" — 1 as r — cxd and it plays no further role in our 
analysis. 
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-4 



Figure 5. ci (1^^ as a function of mass M and cliarge Qi for stable su(2) black 
hole solutions with A = — 10 and varying r^. Only solutions for which the gauge field 
function has no zeros and a;(r/i),ajoo > 1 are plotted. From this surface it is 
evident that ci can be regarded as a single- valued function of M and Qi. 
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Figure 6. Plot of M against charge Qi for su(2) black holes with A = —10 and 
varying r^. There is no evidence that the curves for different intersect, suggesting 
that M and Qi uniquely characterize the black holes. 




( 130|) are satisfieq]]!. We would like to show that the Ci and C2 are determined by M, Qi 
and Q2, where the non-Abelian charges Qi and Q2 are given by ( l48l) : 

= 1 - + ^ul^, = (1 - ^u^loo^ ■ (51) 

II For ease of reference, we call these "potentially stable" black holes. We are not claiming to have 
shown that all black holes satisfying the inequalities (pO| are stable. While the inequalities (|30|) are 
sufficient to guarantee stability in the sphaleronic perturbation sector, they do not guarantee stability 
in the gravitational sector |17j . 
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To produce sensible plots, we first fix M = 10 ± 0.1 (there has to be some range of 
values of M to obtain clear plots). We perform a scan over the black hole solutions, 
numerically integrating the field equations for a grid of values of r/j, Ui{rh), and (X'2(r/i). 
The results are shown in figures [THHl 




Figure 7. Scatter plot of non-Abelian charges Qi and Q2 for potentially stable su(3) 
black holes with A = —3 and mass M = 10 ± 0.1. We perform a scan over the black 
hole solutions, and each data point represents a numerical black hole solution. The 
discreteness of the grid we use can be seen as the values of the charges become more 
negative. We find bands of charge values, although the bands seem to merge and the 
structure becomes less clear as the charges become more negative. 

The values of Qi and Q2 for these black holes are shown in figure [71 We observe 
that the values of the non-Abelian charges Qi and Q2 have a band-like structure, with 
the bands curving round the origin. As the charges become more negative, the structure 
becomes less clear due to the finite grid on which we find our solutions. We found similar 
behaviour for other values of M. 

Scatter plots of ci and C2 as functions of Qi and Q2 for M = 10 ± 0.1 can be found 
in figure |8l It is difficult to see clearly from figure El but rotating the scatter plots 
on a computer screen indicates that both ci and C2 are single-valued functions of Qi 
and Q2- Again, we find similar behaviour for other values of M. This provides some 
(albeit limited) numerical evidence for the parameters Ci and C2 being determined by 
the non-Abelian charges Qi and Q2- 

As in the su(2) case, further numerical evidence that the mass and non-Abelian 
charges characterize the black holes can be found by plotting mass M as a function of 
the charges. Firstly, in figure |9l we plot M as a function of the effective charge Q f H6|) . 
to see if the quantites M and Q are sufficient to uniquely fix the black hole solution. In 
figure ini we see that the surfaces corresponding to different values of the event horizon 
radius r^ overlap. For example, there are black holes with mass M = 6.1 and effective 
charge Q = 5 for rh = 1, 1.25 and 1.5. We therefore conclude that su(3) black holes are 
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Figure 8. Scatter plot of the constant ci (left) and C2 (right) against the non-Abelian 
charges Qi and Q2 for the same black holes as in figure [T] It appears to be the case 
that ci and C2 are single- valued functions of Qi and Q2- 
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Figure 9. Scatter plot of mass M as a function of effective charge Q (|46|) for su(3) 
black holes with A = —3 and event horizon radius = 1 (blue), 1.25 (green), 1.5 
(red). It can be seen that the surfaces for different rh overlap, indicating that M and 
Q do not uniquely characterize the black hole solutions. 

not uniquely specified by their mass M and effective charge Q. 

On the other hand, if we plot M as a function of the non-Abelian charges Qi and 
Q2 (ISTj) . a different structure can be seen in figure [TOl The surfaces in figure [TO] appear 
to foliate the parameter space, and we found no evidence of the surfaces for different 
rh intersecting. This indicates that the parameters M, Qi and Q2 do indeed uniquely 
characterize the black hole solutions. 

5. Analyticity of M , Qj and cj as functions of horizon parameters 

The numerical analysis of the previous section has indicated that the + 1 parameters 
M, A and Qj, for j = 1, . . . N — 1, completely determine stable black hole solutions. In 
section Ej we will give an analytic argument to show that this assertion holds, at least 
for a large subspace of stable black holes. The essence of the argument is to consider 
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Figure 10. Scatter plot of mass M as a function of non-Abelian charges Qi and Q2 
([51]) for the same black hole solutions as in figure [9] The different coloured surfaces 
are for different values of the event horizon radius r^- There is no evidence that the 
surfaces for different intersect, suggesting that M , Qi and Q2 uniquely characterize 
the black holes. 



the + 1 parameters (A, M, Qi, ■ . . ,Qn-i) as functions of the the + 1 parameters 
{A,rh,ijji{rh), . . . ,uJN_i{rh)) which are known, from [16], to uniquely characterize black 
hole solutions of the field equations. If the map 

(A, r^, uji{rh), . . . , UN-iirh)) (A, M, Qi, . . . , Qn^i) (52) 

is invertible, then we can deduce that (A, M, Qi, . . . , Qat.i) uniquely characterize the 
black hole solutions. A convenient way to show that the map ( 152|) is invertible is to use 
the inverse function theorem, so that if we can show that the Jacobian J, given by 

1 



J 














dM 


dM 




duN~i{rh) 


dQi 




duJi{rh) 


dujN-i{rh) 


OQn-i 


OQn-i 


dui{rh) 


duN-i{rh) 


521) is invertible. However 



(53) 





dM 

drh 
dQi 

drh 



drh 



flS2]) is continuously differentiable. In fact, in this section we will be able to prove a 
stronger condition, namely that the quantities M, Qj and cj (122|) are analytic functions 
of {A,rh,uji{rh), ■ ■ .,UN-i{rh)). 

We begin by writing the field equations f lTOl [T3|) in an alternative form, following 
Proposition 7 in jl6]. New variables A, ipk and are defined by 

Ar^X , ^ r"^ duj , ^, 

- ^ - ^ , ^j, = — — , (54) 



rh 



rh dr 
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and then the field equations take the following form: 

dX dijjj dij 

^TTT = ^/a, z— = -z^„ z— = zf^^ 



4 = ^'fs. 4 = 0, .^ = 0, (55) 
dz dz dz 

where the new independent variable \s z = r^/r, so that z G [0, 1]. The functions /a, f^- 

and fs can be found in [16] (noting that in that reference, is rh times as defined 

in (15^ ): their detailed form is not essential here, their key property being that they are 

analytic functions of z, A, r/j. A, the ipj and the ^j, at least in a neighbourhood of 2; = 0. 

Proposition 7 of [16] then gives local existence of solutions of the differential equations 

flSSjl . in a neighbourhood of z = 0, these solutions being analytic in M, A, r^, uj^oo and 



Cj/rh. 



The equations (155|) are regular only on the interval z G [0, 1). In the proof of 
Proposition 6 in [16j, equations of a similar form to fl55l) . but with different independent 
and dependent variables, are derived and used to show the existence of local solutions 
in a neighbourhood of the event horizon, and analytic in A, rh, and ujj{rh). These 
latter equations are regular on an interval equivalent to 2; G (0, 1]. To show that M, 
ojj^oo and Cj are analytic functions of rh, A and ujj{rh) we therefore have to match the 
solutions which exist in a neighbourhood of the event horizon with those that exist in a 
neighbourhood of infinity. 

To do this, let 5 = 2; — |. Then the equations ( l55l) can be written in terms of this 
new independent variable as 

,§ = i.^fs. 4 = 0, 1^ = 0. (56) 

dz dz dz 



Therefore, by analogy with fl55|) . we have a local existence theorem for solutions of the 
above differential equations, in a neighbourhood of z = (2; = |), and the solutions are 
analytic in rh. A, and the values of A, S, ipj and at 2 = 0. 

Since the equations (156|1 are obtained from ( 155|) by a simple translation of the 
independent variable, the new equations (|56l) are regular at 5 = — |, which corresponds 
to z = 0. Therefore the values of the field variables at z = are analytic in rh, A and 
the values of A, 5*, tpj and at 5 = 0. In particular, the quantities M, Qj and Cj will 
be analytic in rh, A and the values of A, S, ipj and at z = 0. 

Now suppose we have a space of black hole solutions of the original field equations 
(fTOl [T3II which are regular everywhere between the event horizon and infinity. In 
particular, every member of this space of solutions is regular at z = (infinity), z = ^ 
and 2; = 1 (event horizon). From Proposition 6 in [16], these solutions are analytic in rh, 
A, and ujj{rh) in the neighbourhood of the event horizon in which they exist. Therefore, 
the values of the field variables at 2; = ^ (2 = 0) are analytic in rh, A and Uj{rh)- 

We now have the following situation: the values of the A, S, tpj and at z = 
are analytic functions of rh, A and Uj{rh). In addition, the quantities M, Qj and Cj are 
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analytic in r/j, A and the values of A, S, ipj and at z = 0. Therefore we can conclude 
that M, Qj and Cj are analytic in r^, A and the Uj{rh). Furthermore, as a corollary, it 
must be the case that the Jacobian J ( l53i) is also an analytic function of r/^, A and the 
(^j{rh)- 

6. Characterizing su(A^) EYM black holes for large |A| - analytic work 

Our purpose in this section is to argue that the charges Qj psj) . together with the mass 
M and cosmological constant A, are sufficient to uniquely characterize 3u{N) hairy black 
holes, at least for some subset of those black holes which are linearly stable [17j . 

There are a number of difficulties in deriving such a result: firstly, general black 
hole solutions are known only numerically; secondly, while the local existence theorems 
in [I6j near the event horizon and at infinity are valid for any values of the cosmological 
constant, event horizon radius and the other parameters in the theory, the theorems 
proving the existence of regular black hole solutions are valid only for sufficiently large 
|A|. 

In other words, in [16] the existence of black hole solutions is shown for fixed r^ 
and Uj{rh), and sufficiently large |A|, where, as discussed in [16], how large "sufficiently 
large" is will likely depend on the values of the parameters r^ and Uj{rh). In the present 
analysis we wish to study the problem from a different perspective: we want to fix |A| to 
be some (suitably large) value, and vary rn and Uj^rt), and then show that, at least for 
some subset of the black hole solutions thus generated, the black holes can be uniquely 
characterized by their mass M and non-Abelian charges Qj. 

From our above discussion of the generalized "no-hair" conjecture, it is clear that 
we are only interested in characterizing stable hairy black hole solutions. Due to the 
symmetry (fT5l) of the field equations, we may consider only ujj{rh) > without loss of 
generality. We therefore restrict attention to black holes for which the inequalities fl30l) 
are satisfied, both on and outside the event horizon. It is straightforward to show in 
this case that the functions 0Jj{r) are monotonically increasing and that ujj{r) > for 
all r >rh- We also consider only those black holes which are thermodynamically stable 
(see section [23]) . which is equivalent to focussing attention on comparatively large black 
holes. 

6.1. Characterizing black holes by M and Qj 

We begin our analysis with a more careful consideration of the space of black hole 
solutions for large |A|. There are two length scales in this problem: the event horizon 
radius r^ and i = a/— 3/A. We wish to study black hole solutions when £ is small, 
but allowing r/j to vary. To this end, it is helpful to introduce a dimensionless radial 
co-ordinate x = r/rh, so that x G [1, oo) irrespective of the value of the event horizon 
radius r/j. In terms of x the field equations (fTOt [T3!) take the form: 
dm „. on 
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2 d'^Uj 

dx'^ 



+ 



2m — 2r?x'^r 



+ 



2r\x^ 



dujj 
dx 
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(57) 



where m{x) = 
take the form 



m[r 



)/rh, the quantity Wj is given by ( IT2|) and equations (fTT| [T4|) now 



TV 



4^,4^4 



2 2 



N-i + 2jy 



1 y. f dujj 
r? ^ V dx 



G 



At the event horizon, x = 1, we have 



mil) 



- 1 + — 

-I ^ £2 



(5^ 



(59) 



which becomes large as £ — )■ 0. We define a further new variable fh{x) by 

m(x) = mi + m(x) (60) 

where mi = m(l) (l59l) . In ^16j, it is shown that ^ (m£^) — )■ as £ — for fixed r/^, 
Uj{rh), but this does not necessarily mean that ^ tends to zero in this limit. Indeed, 
we shall see below that this is not the case, as is borne out by Figures 5 and 6 in [T6] . 
We now write the first and third field equations f l57|) in the form 



^2dm 1 
dx r? 



2mii 



2mf 



n N-1 



X 



X 



I 2 2 
+ V 



AT 



+ ifvE|(-l--l-.-^-i + 2^) 



2mi£^ 2mf 



X 



X 



I 2 2 
+ V 



+ [2mi^2 ^ 2m£2 _ 2rlx^fpe 



d'^Uj 
dx"^ 



(61) 



The analysis leading to Proposition 11 in [16] essentially involves ignoring the terms 
mP' and i'^WjUj in the right-hand-sides of the above equations. It is clear that for any 
fixed values of r^ and Uj{rh), we may choose £ sufficiently small that these terms are 
negligible compared with the others in the equations. However, here we wish to vary 
r/i and Uj{rh), and therefore we need to carefully examine the magnitudes of all the 
quantities in the above equations for small i. It is clear from equations (16T]) that there 
are subtleties if rt is also small. For this reason, we consider only those black holes for 
which r/i ^ i, as these will be thermodynamically stable (see section [275]) . 

For a regular event horizon at x = 1, we can vary rh and ujj{rh) subject to the 
single constraint ( 12T|) . which can be written as 

N 



EN' 



rh 



N-l + 2jy 



<2rlf + 6rl 



(62) 
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< 2rlf 



+ 6ri 



(63) 



Therefore, for each j, we have 

f (a;|(r,)-u;|„,(r,)-iV-l + 2j) 
and, as a result, it must be the case that 

i [Lo]{r,) - J (iV - j)] < J [2rlf + 6r^] ^ . (64) 

Our numerical results in section 12.41 indicate that we do not have regular black hole 
solutions for values of (jJj{rh) close to the boundary of the region defined by (|62l) . and 
therefore we do not need to consider all ujj{rh) such that fl62|) is satisfied. At the same 
time, we find that the region of parameter space for which we have regular black hole 
solutions which are uniquely determined by M and Qj seems to grow as £ decreases. 
Therefore we would like to consider a region of the ujj{rh) parameter space which, for 
small i, is smaller than the region defined by the inequality ( 162|) , but which nonetheless 
grows as i decreases. 

To this end, define new functions qj{x) by 

f[u;^(r)-j{N-j)Y = f-q^{x), (65) 

where a > is a presently unknown constant, and is the same for all j. We further 
assume that qj{x) is order one for small i. Setting a = 1 would correspond to an 
upper bound on Uj which is fixed independent of i, whilst setting cr = corresponds 
to considering the whole of the region of parameter space satisfying ( l62l) . We therefore 
anticipate that < cr < 1 will be relevant for our analysis. In this case, there is an 
upper bound on uj of the order of which grows as i decreases. 

In the su(2) case, it is shown in |12] that, for fixed and w(r/i), i^^u'^r) — )■ 
as £ — )■ 0. This suggests the following definition of further new functions rij{x), again 
expected to be order one for small i, such that 



dx 



7]j[X) 



(66) 



for some k > 0. 

Let us now examine whether it is possible to find suitable values of a and k so as 
to obtain consistent approximate solutions of the field equations (|6T1) when i is small. 
First write the equations ( 16T1) in terms of the qj and rjj-. 



,dm 
dx 



2mi£^ 2m£2 



X 



X 



I 2 2 
+ V 



X] 



()2a 



+ 



N 



\qj{x) - qj^x{x)\' 



(67) 



2/)K + l 



x't 



I 2 2 

+ r.x 

ff2 , n^/)2 o„2^3/)2^ 



+ e^+^ [2mi£2 + 2mf - 2rlx^ 



drjj 
dx 

+ 2r2x3] r]j{x) 



+ - 2gi(x) + qj-iix)] uj^ix). 



(68) 
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Examining first (1681) . and bearing in mind the upper bound P'^~^q'j on u}j{x), for non- 
trivial solutions we require that the first two terms are of the same order in i as the last 

term. This means that 

3cr 1 , , 

« = y - 2- m 

Requiring that n > then implies that a > 1/3, which is consistent with our 
assumptions. With this value of k, we then have 

2/? + 2 = 3a + 1 > 2a (70) 

for all a > 0, which means that the first line of (1671) is of subleading order in i compared 
with the second. Futhermore, differentiating fl65|) gives 

^ ~ 2g,,/^+- (71) 

for Uj ~ and k given by fl^ . Therefore we may regard the functions qj{x) as 

constant to leading order in i. Integrating (1671) then gives, to leading order in i, 

N 



fm(a:) = ^(l-i)5^[g,(l)-g,_i(l)f 
1 \ ^ 

1 - - j E - -N-1 + 2jY], (72) 



N 

E 

i=i 

where we have used the initial condition m = at the event horizon x = 1. Note 
that the answer (1721) implies that the terms in (l68l) are indeed small compared 
with i'^rui = 0{1) for small £ and so can be ignored to first order. However, it is not 
necessarily the case that m itself is small. In fact, we have m ~ (9 (£^°"~^), so m will in 
fact diverge as £ — )■ 0, albeit more slowly than mi (IMJl) . 

To leading order, the Yang-Mills equations ( 1681) become 



= x'^ 



2 2 2mi£^l dn 
ry^x 

X 



^ + [2mie + 2rlx'^] r],{x) 



+ \ fe+i(a;) - 2qj{x) + qj-i{x)] gj(x)^ (73) 
where we have ignored the term 

TV 

2rlxH'pe = 2rlxH^^ fe(^) " ^^"i^^)]' ' (^4) 

and used the leading order approximation 

CJ. =£|(— (75) 

To leading order in i, we can treat the functions qj as approximately constant, and in 
this case (173|) can be integrated to give 

'^'^''^ = ' 2rl{x^lx + l) " ^^^'^^^ + ^^-^^^^^ ^^^^ 
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where we have chosen the arbitrary constant of integration to be such that r]{x) is finite 
at the event horizon x = 1. Restoring the original variables, it is straightforward to 
check that the solution (176|) for rjj satisfies ( l20l) at the event horizon. 

Therefore we have obtained a consistent, approximate set of solutions of the field 
equations which are valid for all r/j S> ^ and all u}j{rh) such that 

[oj^(r,)-jiN-j)Y <i"^-' (77) 

for some a G (|, l). 

At this stage we compare our approximate solutions qj ~ constant, rjj fl76p and rh 
(!72|) with numerical solutions to test the validity of our approximations. To illustrate 
the behaviour, in figures [TTHTSl we consider su(3) black holes with r^ = 1, (jJiir^) = 1, 



3 and varying cosmological constant A = —10^, —10^ and —10^. 
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Figure 11. Gauge field functions uJi{r) (left) and aj2(?') (right) for su(3) black holes 
with Tfi — 1, uji{rh) = 1, uj2{rh) = 3 and varying cosmological constant A — —10'*, 
— 10^ and —10^. The gauge field functions approach constants as |A| increases. 



Firstly, in figure [TT] we plot the gauge field functions uji{r) and a;2(r). It is clear 
that, as expected, the gauge field functions approach constants as |A| — )• oo. 

To test the validity of the approximate solution for rjj (176!) . for each value of r 
we divide the numerically-generated value of ^(r) by rjj given by ( |76l) . The results are 
plotted in figure [T2| where we have also divided the answers by a constant (corresponding 
to the powers of i in the definition fl66l) ). so that all curves pass through ±1 at the event 
horizon, to make comparisons easier. Since u}'i{r) < and uj'2{r) > for the particular 
black hole solutions we are considering, we have fixed the additional constant so that 
the curves for the uj'i{r) approximation to pass through —1 at the event horizon, and 
those for the uj2{r) approximation pass through 1 at the event horizon. In figure [T2| we 
see that the curves tend to ±1 for all r as |A| — )■ oo. This means that, in the large |A| 
limit, the approximation rjj ( 176|) becomes increasingly accurate. 

To test the approximation m (!72|) . we take the numerical values of m(r) — m{rh), 
divide by the approximation ( 1721) and, as in figure [I2l we also divide through by a 
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Figure 12. Test of the approximation ((75)) for U![{r) (left) and Lj-'2{r) (right). The 
numerical values of uj'j (r) for the black hole solutions in figure [11] are divided by the 
approximate solution rjj (j76p and a normalization constant, chosen so that the curves 
pass through ±1 at the event horizon. The curves tend to ±1 for all r as |A| increases, 
indicating the validity of the approximation rjj (j76p . 

constant (corresponding to powers of i) so that all the curves pass through 1 as r — > oo. 
We fix the normalization at r — i- oo rather than at the event horizon because, as r — )• r/j, 
both m(r) —m{rh) and m f l72|) vanish and numerical errors become an issue in dividing 
these two quantities. The results are plotted in figure [131 In figure [13] we see that 
the curves tend to 1 for all r as |A| — t- oo, showing that, in the large |A| limit, the 
approximation m ( !72|) becomes increasingly accurate. 




10 100 1000 r 10000 

Figure 13. Test of the approximation (|72l) for m(r) — m{rh). The numerical values of 
m(r) — m{rfi) for the black hole solutions in figure [TT] are divided by the approximate 
solution for m (|72l) . and a normalization constant, chosen so that the curves pass 
through 1 as r — >■ oo. The curves tend to 1 for all r as |A| increases, indicating the 
validity of the approximation for m(r) — m{rh) ()72p . 
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We emphasize an important point here. In our su(3) numerical examples above, 
we fixed r^, uji{rh) and U2{rh) and examined the solutions for varying values of the 
cosmological constant A. However, this is not the purpose of the approximation 
developed in this section. The approximation derived here is for fixed (but very large) 
A, and varying rh (keeping rt^ i), uji{rh) and U2{rh) (subject to qj ( l65l) being of order 
one for small i). We think of the approximation Uj = constant, rjj fl76|) . m fl72|) as being 
the first term in an asymptotic series for the field variables, which is asymptotic for large 
|A|. This approximation is uniformly valid on a region of the parameter space in which 
r/i ^ i, and a (l65l) is fixed between ^ and 1. This does not cover the whole space of black 
hole solutions for this large, but fixed, value of |A|. However, this is not unexpected: 
the original existence theorems [121 US] were only proved for | A | "sufficiently large" and 
fixed rh and uj. It is therefore unsurprising that we are only able to derive analytic 
approximations on some subspace of the set of black hole solutions. However, the space 
of black hole solutions that we are able to describe is large and increases in size as i 
decreases. 

Having derived the approximate solutions Uj = constant, rjj fl76|) . m f l72|) for small 
i, we now examine whether these approximate solutions are determined by their mass M 
and non-Abelian charges Qi. For these approximate solutions, the gauge field functions 
Uj are approximately constant, and therefore the charges ( HHj) are given, to leading order 
in i, in terms of the values of the gauge field functions on the event horizon: 

In addition, the masses of the black holes are given, to leading order in i, by 

■ 3 1 ^ 

^=y-^ + i^E [(-K-^) - -Ur^) -N-1 + 2jY]. (79) 



To show that the Jacobian J (13^ does not vanish in this case, we require the partial 
derivatives 

dM 1 Arl 1 



- — __ = Wk{rh)uJk{rh); 

duk[rh) r-h 

dQ^ =0; 
drh 



^-i^E[(-J('''.)--j-.('''-)-^-i+2,r 



dQi _ s/iij + 1) 2Mrh) 



The fact that the Qj do not depend on r^ means that the Jacobian J ( l53|) is given by 

^=g^^Q. (81) 

where Jq is the Jacobian of the charges Qj in terms of the ujk{rh). We observe that 
1^ > follows immediately from fl2ip . while Jq must be non-zero by virtue of the fact 
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that the transformation (HSj) from the values of the Uk to the charges Qj is invertible. 
Therefore the Jacobian J f lHT]) is non-zero. This means that the parameters A, M and 
Qj uniquely specify the black hole solutions, at least in this approximation. 

The approximation we have used in this section is exact in the limit £ — )■ 0. We have 
therefore shown that the Jacobian J is non-zero in this limit. Using the fact that the 
Jacobian J is analytic in the parameters of the theory, as shown in section [5l it follows 
that J 7^ at least for sufficiently small i. We therefore deduce that at least a subset 
of stable hairy su{N) black holes, for sufficiently large |A|, are uniquely characterized 
by their mass, the (negative) cosmological constant A and a set of global conserved 
non-Abelian charges Qj. 

A further comment is in order. In previous sections, we have emphasized that 
Bizon's modified "no-hair" conjecture applies only to stable black holes. In our 
numerical work in section HI we restricted our attention to "potentially stable" black 
holes, namely those for which the inequalities ( 130|1 . necessary (but not sufficient) for 
stability, are satisfied. However, in this section we have made no reference to these 
inequalities, which, with the definition f l65|) . take the form 



Simply by restricting our attention to those black holes for which ( !82|) hold at the event 
horizon x = 1, we trivially have that stable black holes are uniquely specified by A, M 
and the charges Qj, for sufficiently small i. 

6.2. Characterizing solitons by M and Qj 

In the previous subsection we have shown that 5\x{N) black holes, at least for sufficiently 
large black holes in the presence of a sufficiently large |A|, are uniquely characterized by 
their mass M, the su(A^) charges Qj and A. One remaining issue is whether it is possible 
for soliton solutions to have the same M and Qj as black hole solutions for a particular 
A. In other words, is it possible to confuse black holes and solitons by measuring their 
mass M and non-Abelian charges Qjl 

Our analysis mirrors that in the previous subsection. However, as in |16], we find 
that analytic work with soliton solutions is considerably more complicated than that 
for black hole solutions. The work in this section follows the approach and notation of 
[T6] in dealing with the soliton solutions, and we refer the reader to that paper for more 
details, keeping the presentation in this section brief. 

Our goal is to find an approximation for the soliton solutions which is valid for 
small I. Since there is no event horizon for soliton solutions, we have just one length 
scale, namely and we define a new dimensionless radial variable y hj y = r/i. Our 
first task is to write the field equations (fTOl [T3l) in a form suitable for analysis. Following 
[IE], we write the gauge field functions ujj{r) as 



Qji^) > 2 fe+l(^) + ^i-l(^)]• 



(82) 



^iW = [3 {N - 3)f^ Ujiv) 



(83) 
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and define a vector u = {ui, . . . ,un^i)'^. As in [TB], we rewrite the vector n as a sum 
over eigenvectors of the (A^ — 1) x [N — 1) matrix A whose entries are: 

= [j (N - j)]^ [25,, - - , (84) 

where 6ij is the usual Kronecker delta. The form of u in terms of eigenvectors of A 
reads (cf. (77) in [16]): 

TV 

uiy) = uo + J2f3kiy)y'^'' (85) 



k=2 



where Wg = (1, 1, . . . , 1) and the (3j^ are vector functions which satisfy 

A(3, = k{k~l)(3,. (86) 

Next we define scalar variables Ck{y) by (cf. (81) in [TB]) 

a{y) = (Tlf3,{y), (87) 

where cr^ is the k—th left eigenvector of the matrix A. 

We next define a rescaled metric variable rh{y) as follows: 

My) = (88) 

which satisfies the Einstein equation (fT3|l : 

dm ^ I ~ , 

where G is defined in f |T^ and 

P{y) = 2y^e^pe, (90) 

with pg given by (fTTI) . 

It is shown in [16] that the soliton solutions are determined, in a neighbourhood of 
the origin, by the N parameters A (or equivalently, i) and Cfc(0)l3ll. Unlike the black hole 
solutions, there are no a priori bounds on the values of Cfc(O) for the existence of regular 
solutions. However, numerical analysis [I5] shows that the region of parameter space for 
which soliton solutions in which all the gauge field functions Uj have no zeros expands 
as i decreases. It is argued in [31] that, for su(2) solitons, the size of the parameter 
space expands 3jS £ clS i decreases. With this in mind, following (|65|) . we define new 
variables ak{y) by 

Ck{y) = akiy)n-\ (91) 

where each cxfc is a constant, and we expect that 0"^ < 1 so that the space of soliton 
solutions that we are considering expands as i decreases (a^ = would correspond to 
the results for su(2) solitons [31], but, as with the black holes, we expect to not be 
able to consider the entire solution space, but nonetheless a region of solutions which 
expands as £ decreases). Note that, unlike the black hole case fl65l) . we allow for the 

% As with the black hole solutions, strictly speaking there is an additional parameter, namely 5(0), 
but we do not need to consider this further. 
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possibility of different values of ak for different (k- We will assume that each ak{y) is 
order one for all y E [0, oo) and small i. 

The Yang-Mills equation ( |TOl) now becomes 



= l/> 

+ 



y'^ + 2ky'~'^-^ + kik-l)y'-W 



dy"^ 

2m + 2y^ -P 



y 



dy 

k do^k 
dy 



+ ky'' 



+ 



1 



(92) 



where P is given by f lUU]) . and the vector W is defined as W = (W^i, W2, ■ 
with the Wj given by ( !T2|) . 

To determine the leading order (in i, for small i) behaviour of the Yang-Mills 
equation fl92|) . we first consider the term cr'^W. It is shown in [16] that this takes the 
form 

z 

<tIW = -k{k- 1) yH^^^^^-^a^ + alT^yH\ (93) 

j=fc+i 

for some Z G N. The cr^Tj are rather complicated expressions which involve products 
of up to three of the Cfc- Since each < 1, the leading order behaviour of cr^r^ 
will be from terms involving products of three Cfc? which will be of order ^i-s+o-a+o-s+o-c 
for some a, 6, c. These will be subleading compared to the first term in ( 193|) if 
^ + cTfc — 1 < j — 3 + (Ta + cTfe + (Tc, bearing in mind that j > k. This inequality is 
satisfied if 1 > a-, > | for all j. Therefore, to find the leading order behaviour of the 
Yang-Mills equation fl92|) . we need to keep only the first term in f l93|) . This gives, to 
leading order in 



= /i 



y^ 



yd^ak 



dai 



+ 2ky—^ + k{k-l)ak 
dy dy 



2m + 2y^ - P 



dak 
dy 



+ ky 



— k {k — 1) ttk- 



(94) 



Further simplification of the equations fl89| [94]) requires analysis of the quantities 
G and P. First we write the vector functions in terms of the scalar variables (k fISTI) 
(with no summation implied): 

Pkiy) = VkUy) = vd'^^-^auiy), (95) 

where the = (f Vk,2, ■ ■ ■ , Vk,N-i) are right-eigenvectors of the matrix A (JHlj). Using 
the expressions fl83l [85]) , we find 

N 

l + Y^Vkyi'^-'-'-'akiy) (96) 



k=2 



and therefore the leading order behaviour of G (1141) is 

-I 2 



G = £2-2 S 



G 



y 



,da2 
dy 



+ 2ya2 



(97) 
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where 

N-l 

T.G = Y.j{N-j)vl^. (98) 

i=i 

From ( l90l l96|) . it can be seen that P is a comphcated sum of terms involving products 
of at least two Cfc, or, equivalently, at least two a^. Products of the form afciCt^j in P 
are multiplied by coefficients of order ^'=i+^2+o-fei+o-fc2-4 ^^^^^ therefore the leading order 
term of this form occurs when ki = 2 = k2 and is of order i'^'^'^ . Products of the form 
ttfcittfc2Q^fc3 are multiplied by coefficients of order f'^i+'^z+fes+o-fc^+o-fc^+o-feg-s 
leading compared with the 0{i'^''^) term. Similarly, products of the form ak^ak2(yk30:k4, 
are also sub-leading, and we deduce that the leading order behaviour of P is 

P = 2f"2Spa2^3 + o(£), (99) 

where 

N-l 

Ep = ^ [j {N - j) V2, - U - 1) (iV - J + 1) V2,^^f . (100) 
i=i 

Substituting for G and P in the Einstein equation (l89il . a consistent, non-trivial 
solution exists when 

rh = f'''x{y), (101) 

with x{y) satisfying, to leading order in i, the differential equation 



2da2 
?/ — + 2ya2 



1 2 



+ Spa^y^ (102) 



dy 

The m and P terms in f l94p can therefore be ignored to leading order in i provided 
a2 > 0, giving simplified Yang-Mills equations 

= y{l + y^)^ + 2[k + {k + l)y^]^ + k{k + l) yak. (103) 

The Yang-Mills equations fll03p were also derived in [16] in the limit £ — )■ 0, but here we 
have used a more subtle approximation (by including the ak in (19T1)). Equations f ll03p 
have the following solution regular at the origin |16] : 

My) = 2F1 Q [A; + 1] , |; A; + ^; «fc,o, (104) 

where 2P1 is a hypergeometric function and = afc(O). The function ( 1104p has 
magnitude bounded by |afc,o| and, as ?/ — > 00, it tends monotonically to zero as 0(y~^'^"'"*^) 
for k > 2, and as 0{y~'^) for k = 2. Then ( I102p can be integrated, to give x(2/)) which 
is then a bounded function satisfying the required boundary conditions: 

X(l/) = 0(i/3), y->0; X(Z/) = Xoo + 0(i/-i), y^oo. (105) 

We therefore have a consistent set of solutions to the field equations, valid when i <^ 1 
and providing | < cr/j < 1. 



Characterizing black holes with abundant gauge field hair 



31 



Returning the original variables m(r) and ujj{r), we have m(r) = O {P'^^^^) for all 
r, so that the mass of these soliton solutions is also 0{i'^'^^~^^). On the other hand, the 
charges ( l48l) can take on large values because 

r N 1 



and we have < 1. Therefore, for small i, the soliton solutions all have negligibly small 
mass (this is confirmed by numerical calculations). This is in contrast to the black hole 
solutions considered in the previous subsection, which have non- negligible mass provided 
rfi ^ i- We therefore conclude that the soliton solutions cannot be mistaken for stable 
black hole solutions by measuring the mass and non-Abelian charges. 

7. Conclusions 

In this paper we have revisited stable furry black holes in su{N) EYM with a negative 
cosmological constant A, examining the consequences for the "no-hair" conjecture. 
These black holes are "furry" because they have potentially unbounded amounts of 
stable gauge field hair, corresponding to — 1 gauge degrees of freedom for 5u{N). 
Bizon's [8] reformulation of the "no-hair" conjecture states that, within this fixed matter 
model, stable black holes should be uniquely characterized by their mass and a set of 
global charges. Our purpose in this paper has been to investigate whether this modified 
"no-hair" conjecture holds for 5u{N) EYM. 

The space of black hole solutions of su{N) EYM in anti-de Sitter space is 
extraordinarily rich. As well as purely magnetic, spherically symmetric, black hole 
solutions (the focus of our work in this paper), there are also dyonic spherically 
symmetric black holes with both electric and magnetic charges |Tl]; static axisymmetric 
solitons [32] and black holes [33] and rotating black holes [24], as well as a plethora of 
soliton solutions of each of these classes [21]. Furthermore, there are black holes with 
non-spherical event horizon topology [19]. These solutions are known only numerically, 
and, given the complexity of the field equations, deriving general analytic results for all 
these cases is challenging. For general A^, it is known that for sufficiently large |A|, there 
exist purely magnetic, spherically symmetric, soliton and black hole solutions which are 
stable to spherically symmetric, linear perturbations [H | [T6 | [T7]. 

We do not claim to have made a comprehensive study of all families in this 
zoo of solutions. Instead we have focussed on purely magnetic, static, spherically 
symmetric black holes. Within this simplified, restricted model, we have been able 
to present numerical and analytic evidence that black holes which are stable both 
thermodynamically and under linear, spherically symmetric perturbations, are uniquely 
determined by their mass M and a set of A^ — 1 non-Abelian global charges Qj, 
constructed following [28]. Our analytic argument is based on the first term of an 
asymptotic series for the metric and gauge field functions, valid in the limit as the adS 
radius of curvature £ — 0. Furthermore, this series is only applicable for a subspace of 




(106) 
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the full space of black hole solutions in this limit: in particular, we restricted attention 
to large black holes with event horizon radius r^ ^ i and a (albeit large) subset of the 
parameter space of values of the gauge field functions on the event horizon. 

With these limitations, we have provided evidence that Bizon's modified "no-hair" 
conjecture holds for at least this subclass of furry black holes in su{N) EYM with A < 0. 
Of course, it would be of great interest to extend these results to other families of black 
hole solutions within this matter model, particularly dyonic black holes and rotating 
black holes. However, we leave these problems for future work. 

Acknowledgments 

We thank Robert Bartnik for helpful discussions. EW thanks the Universities of Monash, 
Calgary and the Australian National University Canberra for hospitality while this 
work was completed. The work of BLS is supported by a studentship from EPSRC 
(UK). The work of EW is supported by the Lancaster-Manchester-Sheffield Consortium 
for Fundamental Physics under STFC grant ST/J000418/1, and by the European Co- 
operation in Science and Technology (COST) action MP0905 "Black Holes in a Violent 
Universe" . 

References 

[1] Ruffini R and Wheeler J 1971 Phys. Today 24 30-41 
[2] Chrusciel P T 1994 Contemp. Maths. 170 23-49 

Heusler M 1996 Black Hole Uniqueness Theorems (Cambridge University Press) 

Heusler M 1996 Helv. Phys. Acta 69 501-28 

Mazur P O 1987 Proceedings of the 11th International Conference on General Relativity and 
Gravitation 130-157 
[3] Volkov M S and Gal'tsov D V 1999 Phys. Kept. 319 1-83 
[4] Winstanley E 2009 Lect. Notes Phys. 769 49-87 
[5] Bartnik R and McKinnon J 1988 Phys. Rev. Lett. 61 141-4 
[6] Bizon P 1990 Phys. Rev. Lett. 64 2844-7 
[7] Straumann N and Zhou Z 1990 Phys. Lett. B 237 353-6 
Strauniann N and Zhou Z 1990 Phys. Lett. B 243 33-5 
Gal'tsov D V and Volkov M S 1992 Phys. Lett. A 162 144-8 
Volkov M S and Gal'tsov D V 1995 Phys. Lett. B 341 279-85 
Hod S 2008 Phys. Lett. B 661 175-8 
[8] Bizon P 1994 Acta Phys. Polon. B 5 877-898 
[9] Brodbeck O and Strauniann N 1994 Phys. Lett. B 324 309-314 
Brodbeck O and Strauniann N 1996 J. Math. Phys. 37 1414-1433 
[10] Radu E and Tchrakian D H 2011 larXiv:1111.0418l [gr-qc] 
[11] Bjoraker J and Hosotani Y 2000 Phys. Rev. Lett. 84 1853-6 

Bjoraker J and Hosotani Y 2000 Phys. Rev. D 62 043513 
[12] Winstanley E 1999 Glass. Quantum Grav. 16 1963-78 
[13] Sarbach O and Winstanley E 2001 Glass. Quantum Grav. 18 2125-46 
Winstanley E and Sarbach O 2002 Glass. Quantum Grav. 19 689-724 
[14] Baxter J E, Helbling M and Winstanley E 2008 Phys. Rev. Lett. 100 011301 
[15] Baxter J E, Helbling M and Winstanley E 2007 Phys. Rev. D 76 104017 



Characterizing black holes with abundant gauge field hair 

[16] Baxter J E and Winstanley E 2008 Class. Quantum Grav. 25 245014 
[17] Baxter J E and Winstanley E, work in preparation 

[18] Balasubramanian V and Kraus P 1999 Commun. Math. Phys. 208 413-428 

[19] van der Bij J and Radu E 2002 Phys. Lett. B 536 107-13 

[20] Kunzle H P 1991 Class. Quantum Grav. 8 2283 97 

[21] Gal'tsov D V and Volkov M S 1992 Phys. Lett. B 274 173-8 

Kleihaus B, Kunz J and Sood A 1995 Phys. Lett. B 354 240-6 

Kleihaus B, Kunz J and Sood A 1998 Phys. Lett. B 418 284-93 

Kleihaus B, Kunz J, Sood A and Wirschins M Phys. Rev. D 58 084006 
[22] Bartnik R 1997 J. Math. Phys. 38 3623-38 
[23] Kunzle H P 1994 Commun. Math. Phys. 162 371-97 
[24] Mann R B, Radu E and Tchrakian D H 2006 Phys. Rev. D 74 064015 
[25] Visser M 1992 Phys. Rev. D 46 2445-2451 
[26] Chan H-M and Tsou S T 1980 Phys. Lett. B 95 395 

Creighton J D E and Mann R B 1995 Phys. Rev. D 52 4569 

Goddard P, Nuyts J and Olive D 1977 NucL Phys. B 125 1 

Kleihaus B, Kunz J, Sood A and Wirschins M 2002 Phys. Rev. D 65 061502 

Sudarsky D and Wald R M 1992 Phys. Rev. D 46 1453 

Tafcl J and Trautman A 1983 J. Math. Phys. 24 1087 
[27] Brandt R A and Neri F 1981 NucL Phys. B 186 84 
[28] Chrusciel P T and Kondraki W 1987 Phys. Rev. D 36 1874-1881 
[29] Lai C, Oh C H and Soo C P 1987 Phys. Rev. D 36 2532 
[30] Corichi A and Sudarsky D 2000 Phys. Rev. D 61 101501 

Kleihaus B, Kunz J and Navarro-Lerida F 2002 Phys. Rev. D 66 104001 
[31] Breitenlohner P, Maison D and LavrelashviU G V 2004 Class. Quantum Grav. 21 1667-84 
[32] Radu E 2002 Phys. Rev. D 65 044005 
[33] Radu E and Winstanley E 2004 Phys. Rev. D 70 084023 



